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I.  INTRODUCTION 


During  the  past  few  decades,  weather  forecasting  accuracy  has  increased  largely  due 
to  the  improvements  in  numerical  modeling  of  the  atmosphere.  Forecasts  of  1  -  3  days  can 
be  made  with  a  good  deal  of  certainty,  allowing  for  necessary  preparations  to  be  made  for 
oceanic  or  continental  winter  storms.  However,  rapidly  changing  events  such  as  explosive 
cyclogenesis  are  still  often  missed  by  the  numerical  models  and  the  human  forecaster. 
This  remains  a  potentially  hazardous  situation  since  people  are  often  caught  unprepared 
with  little  or  no  warning. 

While  the  fundamentals  of  cyclogenesis  are  generally  understood  and  documented,  the 
factors  which  combine  to  create  a  rapidly  deepening  cyclone  are  still  unclear.  Opertitional 
forecast  models  used  for  1  to  3  d:iy  forecasts  have  consi.stently  tended  to  under-forecast  the 
deepening  rate  of  these  storms  (Sanders  and  Gyakum  1980),  although  more  recent  studies 
have  shown  an  improvement  in  this  area  (Sanders  1987).  A  number  of  reasons  have  been 
suggested  to  account  for  this  failure  including  spatial  resolution,  inititd  conditions,  and 
subgrid-scale  parameterizations. 

A  large  number  of  studies  using  numerical  model  simulations  of  explosive  cyclogene¬ 
sis  have  been  conducted.  In  general,  these  studies  fall  into  two  categories.  The  first  group 
includes  near  perfect  model  simuUitions  allow'ing  for  diagnostic  studies.  These  studies 
allow  researchers  to  get  a  more  complete  and  dynamically  consistent  picture  ol  the  de¬ 
veloping  cyclone,  even  though  actual  data  measurements  are  sparse  temporally  and  spatial¬ 
ly  (e.g..  Whitaker  et  al.  1988;  Uccellmi  et  al.  1987).  The  second  kind  ol  study  involves 
investigations  of  model  failures  to  match  the  ob.served  deepening  rate  of  an  extratropical 
cyclone  (e.g..  Anthes  et  al.  1983).  As  numerical  models  usually  under-predict  the  ob¬ 
.served  deepening  rate,  various  parameterizations  and  resolutions  are  varied  in  an  attemjii 
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to  improve  the  model  forecast,  witii  varying  degree.-,  of  success.  This  thesis,  however,  will 
investigate  a  failure  of  a  numerical  model  which  forecast  overdevelopment  of  an  explo¬ 
sively  deepening  extratropical  cyclone. 

Recently,  large  scale  experiments  have  been  conducted  in  order  to  study  the  explosive 
deepening  of  extratropical  cyclones.  The  most  recent  of  the.se  was  the  Experiment  on 
Rapidly  Intensifying  Cyclones  over  the  Atlantic  (ERICA),  conducted  from  1  December 
1988  to  1  March  1989  (Hadlock  1988;  Hartnett  et  al.  1989).  Its  purpose  was  to  provide  a 
large  data  sampling  of  .such  cyclones  over  the  open  waters  of  the  western  North  Atlantic 
Ocean.  Aircraft,  buoys,  increased  rawinsondes,  and  radar  data  supplemented  the  usual 
satellite  imagery  during  Intense  Observation  Periods  (lOP's). 

The  second  storm  of  the  experiment,  IOP'2,  devek^ped  on  l.>-14  December  1988.  The 
Naval  Operational  Regional  Atmospheric  Prediction  System  (NORAPS)  forecast  initial¬ 
ized  ill  12Z  on  12  Dec  88  missed  the  time  and  rate  (d'  rtipid  deepening  of  the  IOP-2  cy¬ 
clone,  although  its  48  h  forecast  of  centrtil  pressure  was  only  6  mb  off.  Using  additional 
buoy  and  drop.scnde  data  obtained  during  ERICA,  this  study  will  investigate  the  successful 
and  unsuccessful  as[)ects  of  the  model's  forecasts. 

A  literature  review'  is  presented  in  section  2  focusing  on  the  hitent  heating  aspect  of 
rtipid  cyclogenesis.  Section  3  summarizes  the  synoptic  events  of  lOP-2  and  briefly  de¬ 
scribes  the  operational  numerical  forectists  (d  this  cyclone.  A  brief  description  of  the 
NORAPS  model  is  given  in  sectiim  4.  The  numcrictil  simultitions  mtide  using  the 
NORAPS  model  are  documented  in  section  followed  by  a  diagnostic  study  of  these 
simulations  in  section  6. 
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11.  BACKGROUND 


Numerous  studies  have  been  conducted  to  investigate  the  causes  of  explosive  cyclo¬ 
genesis  and  the  reasons  for  poor  forecasts  of  these  events  by  numerical  models.  Anthes  et 
al.  (1983)  described  model  simulations  of  the  QE  II  storm.  Their  study  concluded  that 
lower  tropospheric  baroclinic  instability  was  the  major  cause  of  growth  for  that  cyclone. 
This  result  was  confirmed  by  composite  rawinsonde  data  studies  from  weather  ships 
(Rogers  and  Bosart  1986).  However,  Uccellini  (1986)  has  show  n  that  upper-tropospheric 
conditions  w'ere  also  important  in  the  OE  II  storm. 

Anthes  et  al.  (1983)  state  that  latent  heat  release  appeared  to  be  mi^re  important  in  the 
latter  stages  of  development.  Numerical  simulations  with  and  without  latent  heating 
show'ed  little  difference  in  the  12  h  forecast,  indicating  the  lack  of  importance  in  latent 
heating  for  the  early  period  of  development  of  that  case.  While  the  incretised  latent  betiting 
improved  the  intensity  forectist.  the  predicted  motion  of  the  sti'rm  w  tis  slow  ed  due  to  this 
heating,  resuniiig  lu  iaiger  posiinin  eiioi.>.  Recem  .siudic:^  al.so  li.ivc  .>hov.  n  that  the  verlictil 
distribution  of  latent  heating  has  a  significtmt  effect  on  the  deepening  the  storm,  w  ith  ;i 
low  er  height  of  mti.xim.um  Iteating  producing  a  greater  intensification.  Danard  ( 16S5) 
postulates  that  this  is  due  to  a  lower  level  of  non-uiveigeri^c  aiid  ii.C!et...od  ;,.v.  le'.e! 
convergence  and  vorticity  production.  The  heating  of  the  lower  troposphere  also  causes 
destabilization,  resulting  in  ;i  more  rtipid  intensification. 

Danard  (1985)  studied  two  explosive  deepening  events  in  the  eastern  North  Pacitic 
using  an  8-level  large  scale  (190  km  grid  spacing)  primitive  equation  model.  He  too 
crmcluded  that  both  latent  and  sensible  heating  play  a  major  role  in  cyclone  intensification. 
Three  variations  of  the  cumulus  parameterization  scheme  developed  by  Kuo  ( 1  074) 


indicated  that  the  maximum  amount  of  convective  precipitation  was  directly  propt'rtional 
to  the  amount  of  deepening,  although  the  differences  w'ere  not  large  (5  mh  at  most). 

Leslie  et  al.  ( 19iS7)  conducted  model  simulations  of  an  explosively  deepening  cyclone 
off  the  east  coast  of  Australia.  Their  findings  on  the  role  of  precipitation  agreed  with 
Anthes  et  al.  (198.^).  For  their  ca.se  study  it  appears  that  cumulus  convectitm  is  not  impi'r- 
tant  during  the  early  sUiges  of  the  cyclone  formtition,  but  plays  a  decisive  role  in  the  rapid 
intensification  stage.  They  also  tested  the  importance  (tf  other  small  scale  processes  such 
as  effects  of  topographv,  surface  fluxes,  and  sea-surface  temperature  analyses.  The\ 
concluded  that  all  effects  \'were  necessary  for  evclone  development  and  that  none  could  be 
omitted  without  degrading  the  lorecast  sigmfictmtly.  Even  with  2.^  km  hori/iuital  lesolu- 
tion,  they  tailed  t('  achiece  the  (ibserved  deepening  rate. 

Mullen  and  Baumhefner  ( 1988)  pointed  out  thtit  case  studies  <'1  indi\  idual  storms  can 
be  misleading  ;md  thtit  the  results  ctinnot  always  be  genertilized.  'I'heir  studs  used  tin  1  1 
w  estern  Niartn  ILicitic  Ocetin  storm  ensemble  from  the  NCAR  Community  Climate  .Model. 
From  ti  l-^(l-dtiy  simulation  ot  perpetual  Jtmutiry  conditions,  their  results  indicate  thtit  tottil 
ditibtitic  hetiting  tind  btiroclinic  dyntimics  contribute  equtilly  to  storm  deselopment.  .Addi- 
tiontilly.  sensible  and  Ititent  hetiting  tire  each  responsible  for  approximtitels  one-halt  I't  tlic 
deepening  due  to  diabtitic  processes.  They  also  found  thtit  the  deepening  rate  u.is  insensi¬ 
tive  to  the  extict  specdictition  of  precipiititio.i  pt.rtin.eterization  although  the_\  camioned 
nt..  to  rule  out  other  schcme.s. 

Robertson  and  Smith  ( 198,’)  conducted  a  diagnostic  energetics  tinalysis  of  extraiiopi- 
cal  cyclone  development  using  the  Limited  Area  .Mesosctile  Prediction  Ssstem  (L.WIPS) 
from  Drexel  Ums  ersity.  Two  versions  of  the  model,  mc'ist  and  dry.  uere  run  .m  i\s  o  case 
studies.  Their  results  also  show  that  the  baroclinic  aspects  of  development  are  greatly 
enhanced  bv  both  convective  :md  stable  latent  heat  reletise  in  the  model.  In  the  moist 
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l\)recr;sts.  a  substantial  amount  ot'  eddy  kinetie  energy  is  generated  due  tct  ageostrophic 
fltnv  induced  by  diabatic  heating. 

A  comprehensive  study  on  numerical  prediction  ot  explosiv  e  cyclogenesis  u  as  recent¬ 
ly  conducted  by  Kuo  and  Ltw-Nam  ( IWO)  'I  hey  ran  24  h  simulations  ot  9  nine  explosive 
cyclones  over  the  western  Atlantic  rrcean.  varying  the  grid  resolution,  cumulus  parameteri¬ 
zation.  initial  and  boundary  conditions,  as  well  as  surtace  fluxes.  Their  results  indicate 
that  the  most  crucial  model  component  lor  short  range  (0-24  h)  prediction  of  rapid  cyclo¬ 
genesis  is  the  initial  conditions,  followed  by  the  iKui/iMital  grid  resolution,  precipil.ition 
parameterizatitm.  and  lateral  birundary  condiiions.  w  iili  surtace  fluxes  plav  ing  a  minor 
role.  More  interesting  is  their  fituiing  thal  upriclu  convection  actually  degraded  the  fore¬ 
cast  due  tit  an  incorrect  location  of  heating,  (irid  scale  precipitation  (i.e..  non-conv ectiv e 
or  stable)  associtited  with  the  slantwise  ascent  in  the  warm  front  was  crucitil  for  rapid 
development.  Simulations  using  the  Ar.ik.iw a'Scluibert  ( 10~4)  scheme  produced  superior 
forecasts  to  tlti'se  using  Kuo  ( 19~4).  Hut  this  wa''  thought  to  be  caused  t\v  the  fact  that  tile 
Arakaw a-Schubert  scheme  resulted  in  much  lower  convective  ptecgutation  amounts,  thus 
allowing  mi'ie  ot  tlie  tivailable  moiMure  to  !'e  usevi  hv  tlie  stable  |irecipitation  scheme. 
Thev  also  concluded  that  the  Kuo  scheme  w,ix  too  ea-'V  to  activate. 

Most  recentlv.  twn  studies  w  ere  ciMulucted  I'H  c_v clones  from  the  Cienesis  ot  .Aihmtic 
Low  s  [-.xperiment  ((l.Al.L  )  using  the  N(  )K  ATS  model.  1  lou  et  al.  (  !  ran  the  mod. el  on 
the  Iniensiv  e  Ot’serv  .alion  PeriiHl  ( K  )P)  2  c\  clone  w  ilh  and  w  ithoul  tlie  extr.i  data  obtained 
as  a  result  of  experiment.  Since  the  cveiogenesis  and  develojunent  were  along  the  east 
coast  of  the  I’.S.,  the  additional  dat.i  trom  the  exjieriment  resulted  in  onlv  a  limiteil  im- 
jriovement  on  the  NOl^APS  tinaivses  and  lorecasis  le.g..  1  -  mb  in  >,entral  pressure).  In 
genertil,  NOR.APS  predicted  the  cyclogenesis  event  lairlv  well  in  the  short  term  (24  -  ,>(i  h) 
but  missed  the  strength  tinii  timing  of  a  cold  surue  resulting  in  poiir  forecasts  later  in  the 
period.  Secontl.iry  cyclogenesis  at  .’oi  li  was  eorrectlv  lorecast  bv  the  model.  Wash  el  al. 


(IWO)  studied  the  GALE  IOP9  cyclone.  Using  the  additional  GALE  observations,  the 
NORAPS  optimal  interpolation  analysis  reveals  critical  subsyni^ptic  features  which  were 
important  in  the  development.  However,  NORAPS  incorrectly  deepened  the  wrong  low 
due  to  a  poorly  forecast  250  mb  jet  streak.  Central  pressures  were  under-forecast  by  10-15 
mb  after  36  h.  NORAPS  forecasts  using  the  additional  GALE  data  were  better  on  the 
storm  track  but  still  deficient  on  intensity. 

The  current  study  also  investigates  the  importance  of  latent  heating  in  the  develop¬ 
ment  of  an  explosively  deepening  cyclone.  However,  the  regional  model  used  in  this  study 
operationally  over-predicted  the  deepening  of  the  ERIC.A  lOP-2  cyclone.  The  importance 
of  the  model  initial  conditions  also  will  be  studied.  ERICA  data  (aircraft  dropsondes. 
rawinsondes.  and  buoys)  will  be  used  to  verify  the  model  performance.  The  results  should 
yield  a  better  understtmding  of  numerical  modeling  processes  in  forecasts  of  explosive 
cyclone  de\’elopment. 
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III.  ERICA  IOP-2  SYNOPTIC  SUMMARY 


During  the  development  of  the  ERICA  lOP-2  cyclone,  two  disturbances  of  different 
origins  interacted  over  the  western  North  Atlantic  Ocean  to  produce  a  rapidly  deepening 
surface  low  (Chalfant  1989).  Figure  I  presents  the  Fleet  Numerical  Oceanography  Cen¬ 
ter’s  (FNOC)  500  mb  analysis  for  1200  UTC  12  December  1988^.  A  large  trough  lies 
over  eastern  Canada  and  extends  along  the  east  coast  (d'  the  United  States.  East  of  the 
Virginia  coastline,  two  air  streams  merge,  one  origina  'ng  from  Northwest  Canada,  the 
other  from  the  southeastern  U.S. 

By  1.5/0000,  short  wave  troughs  are  apparent  in  e;ich  of  the  air  streams  (Fig.  2).  The 
northerly  wave  is  kK’ated  over  Michigan  and  Wisconsin  with  the  southerly  disturbance 
over  the  southetist  U.S.  Note  the  phtise  lag  (d'  the  -.'^0"C  isotherm  compared  to  the  height 
;md  wind  fiel.i  in  the  northern  trough  {(ner  Minnesiitti).  indictiting  the  likelihood  of  tiddi- 
tiontil  barcK’linic  development.  The  FNOC'  surface  analysis  for  this  time  (Fig.  2)  reveals  a 
cold  front  extending  southwestward  from  ;i  I(hv  near  40°N  50'W  to  the  southern  tip  of 
Florida.  Behind  it.  a  high  pressure  system  is  nun  ing  (di  the  New  England  coastline.  The 
V.  inds  associated  with  the  high  arc  advecting  cold  ;iir  o\  er  the  w  arm  Gulf  Stream  waters. 
Note  that  buoy  41001  at  .55°N  7.’'W  is  rep('rting  an  air  tcm[)erature  td'  lO'C  and  a  sea- 
surface  (d'  20.5 "C. 

A  surface  low  begins  to  de\  elop  ah'iig  the  tail  emi  of  the  cidd  front  off  the  Florida 
coast  after  1.7/0000  in  response  to  the  southerly  short  wave.  The  low  slowly  deepens  and 
moves  northeastw  ard  over  the  next  12  h  ;is  seen  in  Fig.  a.  Meanw  hile  :it  500  mb.  the 
mrrtherly  trough  has  continued  to  develop  and  move  to  the  Ohio  Vtilley  (Fig.  5).  Note. 

..  Ilereaf'ler.  dale  and  lime  will  be  elenoleel  in  a  d.!\  .hoLir  I'  lC  lormai,  e  u.  12  1200, 
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however,  that  this  upper-level  wave  is  not  associated  with  any  well  defined  surface  low, 
although  a  surface  trough  does  extend  from  the  Great  Uikes  to  Kentucky  with  some  reports 
of  snow  (Fig.  4).  By  13/1800  (not  shown),  a  second  surface  low  has  formed  east  of  Cape 
Hatteras  in  response  to  the  northeriy  trough. 

The  operational  objective  surface  analyses  without  ERICA  data  do  not  accurately 
depict  this  developing  cyclone  situation.  Thu.s,  subjective  analyses  by  Chalfant  (1989), 
utilizing  satellite  imagery  and  ERICA  aircraft  and  buoy  data,  are  u.sed.  At  14/0000  (Fig. 
6),  three  lows  appear  with  central  pre.ssures  between  996  and  998  mb.  Aloft  (Fig.  7),  the 
northern  wave  has  moved  off  the  East  Coast  of  the  US.  Note  that  the  isotherms  are  now 
m  tre  nearly  in  phase  with  the  height  ccmtcuirs.  indicating  little  or  no  further  baroclinic 
development  aloft  is  expected.  However,  at  the  surface,  the  rapid  intensification  event 
begins  at  this  time. 

The  surface  analysis  at  14/0600  (not  shown)  shows  the  lows  have  merged  and  deep¬ 
ened  approximately  17  mb  over  6  h.  By  141200  (Fig.  8).  the  low  has  plummeted  to  at 
least  968  mb  (buoy  at  37.5°N  63.5'\\'),  a  drop  of  30  mb  in  12  h.  As  w  as  expected  from  the 
14 '0000  500  mb  analysis  of  heights  and  tem|XTatures  (Fig.  7).  the  height  of  the  upper-level 
short  wave  (5340  m  over  Delaware)  has  slunvii  little  change  over  the  same  12  h  period 
(Fig.  9).  although  a  closed  circulati('n  has  develi'petl. 

Operational  numerical  forecasts  df  the  lOP-2  cyclone  varied  in  accuracy.  They  are 
evaluated  in  detail  by  Chalfant  (1980).  None  could  successfully  resolve  the  multiple  low- 
centers  at  14/0000.  The  NMC's  Nested  Grid  MchIcI  (NGM)  48  h  surface  pressure  forecast 
valid  at  14/1200  (ncU  shown)  indicates  only  a  984  mb  low,  an  error  of  16  mb.  However, 
the  position  forecast  is  only  100  km  in  error.  The  Global  Spectral  Model,  also  run  at 
NMC,  produced  a  much  better  forecast  w  ith  central  pressure  of  964  mb  (not  show  n)  and  a 
position  error  equivalent  with  the  NGM.  The  operational  NORAPS  forecast  (Simulation  1 
in  section  5)  correctly  predicted  the  rapid  deepening,  slnwcing  a  962  mb  low  at  48  h.  but  a 
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position  error  of  approximately  350  km.  However,  as  will  be  shown,  the  impressive  fore¬ 
cast  of  rapid  deepening  at  48  h  by  NORAPS  and  the  NMC  spectral  model  is  actually  the 
result  of  poor  forecasts  at  intermediate  times. 
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IV.  MODEL  DESCRIPTION 


NORAPS  is  a  regional  model  run  operationally  by  FNOC,  Monterey,  CA  over  four 
areas:  the  Mediterranean  Sea,  the  Indian  Ocean,  the  North  Atlantic  Ocean,  and  the  western 
North  Pacific  Ocean.  The  purpose  of  NORAPS  is  to  provide  high  resolution  forecasts  of 
up  to  48  h  for  areas  important  to  Navy  operations  (Hodur  1987).  The  NORAPS  model  is 
easily  relocatable  to  any  area  of  the  globe.  The  grid  spacing  is  variable,  both  horizontally 
and  vertical,  but  is  limited  by  the  speed  and  memory  constraints  of  the  computer  system. 
For  simulations  in  the  current  study,  60  km  spacing  is  used  on  a  109  x  82  grid  providing  a 
model  domain  of  approximately  6480  x  4860  km  (Fig.  10).  The  vertical  distribution  of  the 
21  vertical  sigma  levels  is  shown  in  Table  1. 


Table  1.  NORAPS  model  sigma  levels 


Mtxlel  Level 

Siema 

Model  Level 

Siema 

1 

0.010 

12 

0.600 

2 

O.O.Ta 

\? 

0.700 

3 

0.067 

14 

0.787 

4 

0.1 0.a 

15 

0.850 

5 

0.150 

16 

0.890 

6 

0.200 

17 

0.920 

7 

0.250 

18 

0.950 

8 

0..”00 

19 

0.975 

9 

0..TS0 

20 

0.990 

10 

0.412 

21 

0.997 

11 

0.500 

NORAPS  obtains  initial  data  from  radiosonde.  PIBAL,  .aircraft  and  ship  observations. 
In  certain  regions,  satellite  vertical  soundings  and  cloud-track  w  inds  are  used.  Optimum 
Interpolation  (01)  is  used  for  the  analysis  of  these  data  with  the  previous  1 2  h  forecast  used 
;is  the  first  guess  field.  In  order  t('  "spin  up”  the  nKKlel.  four  12  h  update  cycles  are  run 
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before  a  realistic  simulation  is  performed.  To  begin  the  update  cycle,  the  NOGAPS  (Nava! 
Operational  Global  Atmospheric  Prediction  System)  analyses  are  used  as  first  guess  field.s. 
Neither  NORAPS  nor  the  global  model  NOGAPS  perform  a  moisture  analysis.  Thus,  the 
previous  12  h  moisture  forecast  is  used  as  the  analysis. 

Boundary  conditions  for  the  forecast  fields  in  the  NORAPS  model  run  are  obtained 
from  NOGAPS  forecasts.  NORAPS  uses  one-way  interactive  lateral  boundary  conditions, 
which  create  a  "blend  zone"  near  the  edge  of  the  model  grid.  If  the  global  model  provides 
a  poor  representation  of  conditions  on  the  NORAPS  boundary,  the  forecast  in  the  iruerior 
of  the  domain  will  eventually  be  contaminated.  The  length  of  time  it  takes  for  contamina¬ 
tion  to  occur  depends  largely  on  the  speed  at  which  the  wind  field  advects  the  boundary 
conditions  across  the  domain.  The  longer  the  forecast,  or  the  smaller  the  model  domain, 
the  more  likely  it  is  that  poor  boundary  conditions  may  adversely  affect  the  NORAPS 
forecast.  For  the  present  study,  the  model  domain  was  situated  so  that  its  left  boundary 
was  far  upstream  of  the  incipient  cyclone  in  carder  to  minimize  the  effects  of  the  boundary 
conditions  (Fig.  10). 

Topography  in  NORAPS  is  obtained  from  the  Ntivy  ten  minute  global  data  base.  A 
terrain  enveloping  procedure  developed  by  Wallace  et  al.  (1983)  is  used  to  more  effective¬ 
ly  represent  terrain  effects  in  the  model.  A  land-seti  table  with  ten  minute  resolution  is 
used  by  the  model  to  define  land  and  sea  grid  points.  Ov  cr  water,  the  albedo  is  set  to  0.09, 
the  ground  wetness  is  1.0,  and  the  surftice  roughness  is  computed.  Over  ice.  these  values 
are  0.6,  0.0,  and  0.0002435  m.  respectively.  Monthly  climatol(igic;d  values  are  used  for 
the  land  grid  points,  as  w'ell  as  the  sea-surface  temperature. 

The  model  uses  an  Arakawa  Type-C  staggered  grid  (Arakawa  .and  Ltimb  1 977).  Time 
integration  is  performed  using  a  split-explicit  scheme  developed  by  Madala  ( 1982).  In  this 
method,  the  solution  is  adjusted  lor  the  linear  terms  which  govern  the  fastest  gravity 
modes,  permitting  use  of  a  relatively  large  model  time  step.  All  :id\ection  terms  are 


second  order.  Fourth-order  diffusion  is  included  in  the  interior  with  second  order  diffusion 
on  the  first  row  and  column  of  the  grid.  High  frequency  time  oscillations  are  controlled  by 
applying  a  Robert  (1966)  time  filter  with  the  smoothing  coefficient  set  to  0.15. 

The  parameterization  of  surface  fluxes  in  the  planetary  boundary  layer  (PBL)  are 
computed  using  the  similarity  theory  (Louis  et  al.  1982),  where  matching  solutions  are 
obtained  for  the  constant  flux  layer  and  the  Ekman  layer.  The  scheme  is  first  order  K- 
theory  where  K  is  a  function  of  the  bulk  Richardson  number. 

Since  the  lowest  model  sigma  level  above  the  surface  is  equal  to  0.9975,  a  predictive 
equation  for  ground  temperature  given  by  Blackadar  (1977)  is  used  for  land  grid  points. 
The  radiation  parameterization  follows  the  methods  of  Sasamori  (1968)  for  long  wave 
radiation  and  Arakaw'a  (1971)  for  short  wave  radiation.  The  effects  of  clouds  are  incorpo¬ 
rated  into  both  radiation  parameterizations,  w’hile  short  wave  radiation  includes  a  diurnal 
cycle. 

Precipitation  in  NORAPS  can  occur  from  convective  processes  (sub-grid  scale)  or 
from  large-scale  condensation.  Convection  is  parameterized  using  a  Kuo-type  scheme 
(Kuo  1965;  Kuo  1974),  in  which  precipitation  is  linked  to  low'-level  moisture  convergence. 
With  each  time  step,  the  model  determines  whether  a  grid  point  meets  the  neces.sary  condi¬ 
tions  for  convection.  First,  there  must  be  moisture  convergence  in  the  PBL.  Mathemati¬ 
cally, 

=  [P.lg)  jv  q\dp  >  0 

where  is  the  moisture  convergence,  P^.  is  the  surface  pressure,  g  is  gravity,  q  is  the 
specific  humidity,  and  the  integral  is  performed  over  a  depth  of  the  atmosphere. 

Secondly,  the  value  of  above  the  LCL  must  be  less  than  that  at  the  LCL  (i.e.,  a 
conditionally  unstable  layer  must  exist  above  the  LCL).  The  depth  of  this  conditional 
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instability  must  be  greater  than  300  mb.  Finally,  the  relative  humidity  at  the  originating 
level  must  be  greater  than  80%. 


If  the  four  conditions  are  met,  the  model  transports  heat  and  moisture  vertically  to 
each  level  (subscript  k)  using  the  equations  given  by  Kuo  (1965), 


_ Mt  {Tc  -  T), _ 

j  f  Cp  (Tc  -  T)  dp  +  jfl  (gc  -  q)  dp 


^  M  (gc  -  _ 

^fC,(T,-T)dp  +  jfL(gc-g)dp  ' 

Here,  T  and  q  are  temperature  and  specific  humidity  respectively,  where  the  c  subscript 
denotes  cloud  and  unsubscripted  denotes  the  environment.  Cp  represep*s  specific  heat  and 
L  is  latent  heat.  Integrations  are  performed  over  the  depth  of  the  cloud  (i.e.,  where  (T^,  -  T) 
>0). 

Kuo  (1974)  modified  his  parameterization  with  the  inclusion  of  the  "b”  factor.  This 
regulated  the  amount  of  heating  verses  moistening  that  occurs  in  the  scheme.  The  fraction 
"b"  is  calculated  following 

100% -i? 

~  100%  -  50%  ’ 

where  R  equals  the  integrated  relative  humidity  over  the  depth  of  the  cloud  and  100%  >  R 
>  50%.  Now  the  equations  are  modified  as  follows: 


(1  -  t)  Mt  (Tc  -  T)k 

J/Cp  (Tc-T)dp 


b  Mt  (gc  -  9)k 
j  f  dp 
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The  role  of  "b”  is  to  allow  for  rapid  moistening  initially  when  R  is  small  (large  b).  As  the 
column  moistens,  R  increases  and  b  decreases.  This  enhances  the  latent  heating  and  de¬ 
creases  the  moistening  rate. 

At  the  time  of  IOP-2  (December,  1^)88),  the  NORAPS  model  operationally  used  the 
Kuo  (1965)  sclicnie.  In  April  of  1989,  the  Kuo  (1974)  scheme  was  incorporated.  Thus, 
the  simulations  in  section  5  reproduce  the  operational  forecast  of  December,  1988,  as  well 
as  the  current  operational  version. 

Stable  (or  nonconvective)  precipitation  is  assumed  to  occur  when  a  grid  point  be¬ 
comes  supersaturated.  Moisture  is  condensed  until  .saturation  is  once  again  reached.  The 
condensed  moisture  is  evaporated  into  the  next  layer  beltnv  until  it  becomes  saturated. 

The  NORAPS  mode!  has  been  used  operationally  by  the  U.S.  Navv  since  1984.  Skill 
scores  show  its  forecast  accuracy  to  be  ctmiparable  with  other  numerical  models  (Hodur 
1987).  In  addition,  the  NORAPS  model  has  been  used  extensively  as  a  research  tool  at  the 
Naval  Oceanographic  and  Atmospheric  Research  Ltihortitory  and  the  Naval  Postgraduate 
School.  Recent  studies  (Liou  et  al.  1990;  Wash  et  al.  1990)  have  utilized  the  model  to 
study  cyclogenesis  during  GALE. 


V.  MODEL  SIMULATIONS 


In  order  to  investigate  the  physical  processes  of  this  rapidly  deepening  cyclone, 
numerical  simulations  of  the  cyclone  were  run  using  the  NORAPS  model.  Forecasts  of  48 
h  were  made  beginning  at  12/1200,  12  h  before  the  initial  development  and  36  h  before  the 
rapid  deepening.  Since  the  intensification  occurred  with  significant  convection  (Chalfant 
1989),  the  convective  parameterization  scheme  of  NORAPS  was  varied  with  all  other 
physics  in  the  model  treated  identically.  The  four  convective  parameterization  schemes 
that  were  run  are:  1)  Kuo  (1965)  cumulus  scheme  (which  was  the  operational  version  in 
Dpfpmhpr  1988):  2)  Kuo  (1974);  3)  Kuo  (1974)  with  vertical  eddy  diffusion;  4)  convec¬ 
tive  and  stable  precipitation  turned  off“. 

A.  SIMULATION  1 

Simulation  1  of  the  IOP-2  event  appears  to  begin  well.  A  weak  low  forms  at  3()°N  off 
the  Florida  coast  by  13/0000  (Fig.  1 1).  But  comparison  with  Fig.  3  shows  that  the  forecast 
low  is  northeast  of  the  verifying  position  by  over  750  km.  This  is  a  large  error  for  a  12  h 
forecast. 

Figures  12  and  13  present  the  (a)  1000  mb  .streamlines  tmd  (b)  500  mb  heights  and 
vorticity  for  the  analysis  (12/1200)  and  12  h  forecast  (13'0000)  respectively.  A  weak  500 
mb  disturbance  is  located  over  Florida  in  the  analysis  as  indicated  by  the  8  s'^  vorticity 
contour  (Fig.  12b).  An  examination  of  actual  rawinsonde  data  (not  shown)  indicates  that 
the  disturbance  did  exist  as  a  result  ot  cyclonic  shear  across  northern  Florida. As  this  dis¬ 
turbance  moves  otf  the  coast,  it  becomes  superimposed  over  the  trailing  cold  front,  as 


2.  Hcrcaflcr,  the  Kuo  (1%5)  and  Kuo  (1074)  schemes  will  he  denoicd  KU()6.s  and  KU()74  rcspeeiivclv. 


shown  by  the  surface  streamlines  (Fig.  12a).  This  system  develops  rapidly  during  the  12  h 
period  as  the  500  mb  absolute  vorticity  doubles  from  8  s'^  to  16  s'^  (Fig.  13b).  However, 
there  has  been  no  corresponding  change  in  the  height  contours  to  account  for  such  a 
change  in  vorticity. 

Note  that  the  southerly  short  wave  mentioned  in  section  4  is  still  over  the  southeastern 
U.S.,  as  indicated  by  the  vorticity  maxima  over  Georgia  in  Fig.  13b.  There  is  very  little 
indication  of  the  existence  of  this  weak  shortwave  over  Florida  in  the  actual  500  mb  analy¬ 
sis  for  13/0000  (Fig.  2),  although  by  13/1200  (Fig.  5)  it  appears  near  30°N  65°W,  but 
apparently  no  '■■irface  development  took  place.  In  other  words,  NORAPS  incorrectly 
developed  a  surface  low  well  to  the  east  of  the  southerly  wave  in  the  lOP-2  event.  The 
existence  of  this  low  at  the  surface  is  not  verified  by  observational  data,  although  the  weak 
500  mb  wave  did,  in  fact,  exist.  Comparing  Figs.  13;i  and  13b,  it  is  noted  that  the  forecast 
1000  mb  circulation  is  southwest  of  the  vorticity  maxima  at  500  mb.  Thus,  the  low  is 
actually  lying  under  negative  vorticity  advection.  Closer  inspection  reveals  that  the  low- 
level  convergence  zone  northeast  of  the  1000  mb  low  is  directly  under  the  500  mb  vorticity 
maxima. 

Figure  14  shows  the  12  h  cumulative  convective  precipitation  forecast  by  the  model. 
Note  that  the  precipitation  maxima  lies  directly  over  the  1000  mb  convergence  zone  north¬ 
east  of  the  low  (Fig.  13a).  This  is  not  .surprising,  since  hnv-lev  el  mcristure  ccmvergence  is 
one  of  the  four  necessary  conditions  for  convection  in  the  Kuo  scheme  (see  Section  3). 
The  convective  precipitation  is  also  ccnncident  with  the  500  mb  vorticity  maxima  (Fig. 
13b).  While  some  surface  observations  dn  indicate  .some  shower  activity  in  this  area,  there 
is  nothing  to  substantiate  the  predicted  2.8  cm  of  convective  rainfall.  Thus,  the  first  12  h 
of  the  model  simulation  can  be  summarized  as  follows;  a  weak  upper-level  short  wave 
trough  moved  over  a  low-level  baroclinic  zone  east  of  Florida,  initiating  cyclogenesis. 
Low-level  convergence  associated  with  the  surface  low  triggered  convection  northeast  (d 
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the  incipient  cyclone.  The  model  convection  produces  an  erroneous  vorticity  maxima  at 
500  mb. 

This  is  even  more  apparent  in  Fig.  15  where  the  18  h  500  mb  heights  and  vorticity, 
and  the  12  to  18  h  convective  precipitation  fields  are  compared.  The  model  continues  to 
successfully  predict  the  real  short  wave  over  the  Georgia  coastline  (Fig.  15a).  But  it  is 
obvious  that  the  height  and  vorticity  pattern  to  the  east  is  anomalous.  The  model  predicts 
12.1  cm  of  convective  rainfall  to  occur  at  a  grid  point  during  the  pa.st  6  h  with  a  secondary 
maxima  of  approximately  7  cm  to  the  east.  Qmiparing  this  to  Fig.  15a,  note  the  colloca¬ 
tion  of  the  12.1  cm  rainfall  maxima  and  the  small  ridging  in  the  5640  m  line  at  500  mb 
likely  due  to  the  diabatic  heating.  Once  again,  the  ccmvective  precipitation  lies  over  the 
convergent  region  north  and  east  of  the  1000  mb  low-  (not  shown). 

As  the  upstream  short  wave  in  Fig.  15a  merges  with  the  developing  system,  rapid 
deepening  takes  place  (7  mb  in  6  h,  17  mb  in  12  h.  and  28  mb  in  18  h).  By  14/0000,  the 
model  has  deepened  the  low  to  980  mb  (Fig.  16).  A  comparison  w'ith  the  subjective  analy¬ 
sis  in  Fig.  6  shows  the  failure  of  the  model  lorecast.  Only  one  low  is  predicted  which  is 
16-18  mb  too  deep  and  too  far  north,  although  a  surface  trough  northwest  of  the  low  does 
hint  at  a  secondary  development  as  does  the  c(>nvective  precipitation  (not  shown).  The 
500  mb  forecast  (not  shown)  still  contains  the  intense  northerly  short  wave  trough  moving 
off  the  east  coast  of  the  US.  How'e\  er.  the  initial  low  has  already  developed  too  much  to 
allow'  the  new  short  wave  to  create  a  new  low. 

Finally,  Fig.  17  shows  the  48  h  NORAPS  forecast  at  the  surface  and  500  mb.  The 
forecast  surface  low  is  only  a  few  millibars  too  deep  (assuming  the  968  mb  buoy  is  a  true 
representation  of  the  central  pressure),  but  the  position  error  is  approximately  500  km  (Fig. 
8).  At  500  mb,  the  low-  is  similarly  too  deep  ;ind  too  far  n<'rth. 

In  summary,  it  appear'  that  the  (operational  NOR.APS  model  predicted  both  500  mb 
synoptic  features  reasonably  w'ell;  the  southerly  siK'rl  w  a\  e  from  Florida,  and  the  norther- 
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ly  wave  from  Gmada.  Unfortunately,  the  model  over-forecast  ihe  development  of  a  weak 
low  w'hich  preceded  the  southerly  short  wave.  By  the  time  the  northerly  short  wave  ar¬ 
rived,  the  initial  low'  was  too  w'ell  developed  to  permit  secondary  development.  Instead, 
the  northerly  wave  could  only  contribute  to  the  deepening  of  the  existing  low. 

Several  studies  have  show'n  that  enhanced  precipitation  occurs  in  an  area  extending 
from  the  northeast  to  the  north  and  west  of  a  surface  low  prior  to  the  onset  of  rapid  deepen¬ 
ing  (Krishnamurti  1968;  Johnson  and  Downey  1976;  Kocin  and  Uccellini  1989).  Kuo  and 
Low -Nam  (1990)  found  that  this  warm  frontal  precipitation  w  as  crucial  to  the  explosive 
development  of  the  cyclone.  The  w  arm  front  predicted  by  the  model,  as  seen  in  the  1000 
mb  wind  field  and  the  convective  precipitation  pattern,  is  in  agreement  with  these  findings 

Examination  of  the  model  output  reveals  that  during  the  development  of  the  surface 
low,  the  low  is  located  either  directly  under  the  500  mb  vorticity  maxima  or  upstream  of  it. 
One  would  expect  the  low  to  he  located  dow  nstream  t'f  the  vorticity  maxima  in  the  region 
of  the  positive  vcirlicity  advection.  In  additimi.  convective  precipitation  charts  indicate  that 
the  convective  nature  of  this  storm  contributed  heav  ily  to  its  development.  While  ship 
reports  and  satellite  imagery  verify  that  convective  activity  was  present  (not  show  n),  it  is 
unlikely  that  the  model  simulation  is  accurate  in  this  respect.  To  investigate  tlie  conv  ectiv  e 
aspect  of  the  model  simulation,  diagnostics  of  the  effects  (9  conv  ection  are  presented  in  the 
next  section. 

B.  SIMULATION  2 

A  know'n  deficiency  of  the  KUOo5  scheme  is  that  it  added  too  much  moisture  aloft 
(Kuo  1974).  In  order  to  correct  for  this,  the  "b"  factor  was  added  w  hich  reduces  the  mois¬ 
tening  w'hen  near  saturation  is  reached  (see  discussion  of  convectiim  pttrameieri/.ition  in 
section  .7).  Simulation  2  used  the  KU074  .scheme  to  see  if  the  change  made  any  difterence 
in  the  lOP-2  c'se. 
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Figure  18  presents  the  6  h  cumulative  convective  precipitation  predicted  by  the  model 
valid  at  13/0600.  Comparing  this  with  Fig.  15b,  it  is  apparent  that  the  KU074  scheme  has 
reduced  the  convective  rainfall  amount  by  more  than  505/.  However,  the  surface  and  500 
mb  forecasts  (not  shown)  at  all  times  were  nearly  identical  to  that  of  Simulation  1(1-2 
mb  differences  in  the  central  sea  level  pressure).  The  reason  for  the  nearly  identical  fore¬ 
cast  is  that  the  stable  precipitation  was  greater  in  Simulation  2.  as  shown  by  Table  2.  In 
other  words,  since  the  Knn74  scheme  caused  less  convective  precipitation,  there  was 
more  moisture  available  for  the  non-convecti\ e  scheme.  As  Table  2  indicates,  the  total 
amount  of  precipitativin  i.n  Simulation  2  is  nearly  equal  that  (d  Simulation  1.  Therefore,  the 
amount  of  latent  heat  released  m  the  two  simulations  is  nearly  equ.al.  resulting  in  compara¬ 
ble  forecasts.  Danard  (  HhSs)  and  Luni  and  Fisberry  ( 16, s'/)  trbserved  a  direct  ct’rrehition 
between  the  maximum  timount  of  convective  precijutation  and  the  pretiicted  deepening 
r;ite.  Kuo  tmd  l.ow-N'am  (1660)  found  that  the  amount  of  stable  precipittition  wtis  key  to 
determining  the  .amount  of  deepening.  In  the  current  study  the  tiual  precipitation  was 
ceirrehited  to  the  deepening  rtite. 

Ttible  2.  Precipittition  (cm)  predicted  by  the  NOR.-M’S  model  for  Simulations  1.  2.  3.  tiiid 
5.  The  timoLints  are  ti  4n  h  sum  ot  till  grid  points  uith  precipitation  ctuised  b_\  the  IOP-2 
cyclones.  .Also  included  is  the  predicted  centrtd  pressLire  (mb)  at  4S  1,. 
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C.  SIMULATION  3 


The  inclusion  of  vertical  eddy  diffusion  in  the  convection  scheme  was  tested  in 
Simulation  3.  The  effect  of  this  diffusion  is  to  mix  the  heat  and  moisture  provided  by  the 
cumulus  parameterization  with  the  surrounding  environment.  This  addition  was  expected 
to  damp  the  effects  of  convection  predicted  by  the  model.  As  seen  in  Table  2,  the  convec¬ 
tive  precipitation  was  further  reduced  by  the  eddy  diffusion.  How'ever,  once  again  the 
stable  precipitation  was  greater  as  a  result.  Thus.  Simulation  3  showed  little  difference 
from  the  previous  two  simulations. 

D.  SIMULATION  4 

.A  fourth  simuhition  was  conducted  in  which  no  precipittition  v.  ;is  allowed.  The  goal 
of  this  simuhition  was  to  ensure  thtit  the  convective  precipitation  scheme  was  truly  respon¬ 
sible  for  the  pocir  forecasts  obserxed  in  the  first  three  simulations.  It  was  unrealistic  since 
no  con\’ecti\  e  or  stable  precipitation  was  allmxed.  In  other  words,  when  a  grid  point 
became  super-saturtited,  the  model  did  nothing  to  correct  this  condition.  Any  development 
that  takes  place  is  due  only  to  adiabtitiv  •■'aroclinic  dynamics,  since  there  is  no  latent  heat 
release. 

Figure  U)  shows  the  (a)  1()(K)  mb  winds  and  (b)  mb  heights  :ind  \orticity  tor  the 
12  h  forecast.  There  tire  mtirked  differences  between  this  simulatii'n  and  Simulation  1 
(Fig.  13).  The  initial  short  wave  at  5(H)  mb  over  f'lorida  does  not  develop  in  Simulation  4 
(note  that  the  enclosed  area  of -S  s'*  vorticity  netir  70  W  in  Fig.  lOh  is  a  minimum). 
Moreover,  no  devehipment  from  this  system  lakes  place  at  1000  mb  (Fig.  10a).  This 
indictites  that  even  the  small  amount  of  convective  precipitation  (i.e.  latent  hetit  release) 
generated  in  Simulation  1  (Fig.  14)  is  necessary  for  the  development  of  the  initial  short 
wave.  Figure  10a  shows  the  beginning  id  the  low  -le\  el  circulation  along  the  Florida  coast 
in  response  to  the  upper-level  short  w  ave  over  (ieorgia.  Thus,  the  mirdel  has  correctly 


ignored  the  initial  w  eak  disturbance  and  developed  a  surface  cyclone  as  a  result  of  the 
second,  stronger  short  w'ave. 

How-ever,  the  model  forecasts  in  Simulation  4  are  in  error  by  not  forecasting  the 
vigorous  development.  The  predicted  surface  cyclone  is  too  shallow  and  too  slow'.  Figure 
20  shows  a  1001  mb  low  at  14/0000  with  an  inverted  trough  extending  along  the  US  coast¬ 
line.  Comparison  with  Fig.  6  demon.strales  that  p’-ecipitation  was  important,  even  at  this 
early  stage.  The  initial  low'  is  3-5  mb  too  shallow  and  475  km  too  far  southwest.  As  one 
might  expect,  the  rapid  deepening  of  the  surface  cyclone  is  completely  missed  by  the 
model  in  Simulation  4.  The  central  pressure  at  14/1200  is  underforecast  by  23  mb  (not 
show'n). 

E.  SU.MMARY  OF  SIMULATIONS  I  -  4 

Figure  21  summarizes  the  four  numeucal  simulations.  It  is  seen  that  the  model  never 
exhibits  a  "rapidly  deepening"  phase,  but  rtither  developed  the  cyclone  at  a  constant  rate, 
although  predicted  deepening  rales  betw  een  13/0600  and  14/0600  (40  mb  in  24  h)  are 
certainly  impressis  e.  In  contrast,  the  observed  cyclone  exhibited  the  rtipid  deepening  in  a 
6  h  period'"’.  The  fact  that,  at  14/0000.  the  first  three  simulations  were  »  15  mb  too  deep 
and  Simulation  4  w  as  only  5  mb  kh'  shallow,  seems  to  indicate  that  the  role  of  convection 
was  rather  small  up  to  this  point.  Still,  the  errors  early  in  Simulation  4  show  that  precipita¬ 
tion  did  play  a  role  in  the  early  stages  of  develcipment.  Oh\  iously.  though,  the  importance 
of  latent  heat  release  in  the  explosively  deepening  phase  cannot  be  ignored.  To  investigate 
this  more  closelv,  several  diagnostic  studies  were  performed  on  the  model  output  and  are 
discussed  in  section  6. 


'ITicrc  is  high  confidenec  in  the  ;ieiu;il  eeninil  pressure  ir;iee  due  to  the  airenit'i  and  buoy  me.'Lsurrrienis 
available. 


F.  OTHER  OPERATIONAL  FORECASTS 


As  mentioned  in  section  3,  the  NGM  and  Spectral  models  run  at  NMC  produced  very 
good  operational  forecasts  in  terms  of  position  error  (100  km  at  48  h).  However,  the  NGM 
was  16  mb  too  shallow  in  its  central  pressure  forecast.  The  Spectral  was  much  better  (4 
mb  too  deep).  But  it  too  was  similar  to  the  NORAPS  model  in  that  it  gave  a  very  poor  36 
h  forecast  (Fig.  22).  As  with  NORAPS,  the  Spectral  model  failed  to  resolve  the  multiple 
centers,  and  the  central  pressure  of  the  low  is  10  mb  too  deep.  The  500  mb  vorticity  fore¬ 
cast  is  not  as  spurious  as  NORAPS,  but  since  a  smoothing  procedure  is  applied  to  these 
operational  charts,  direct  comparison  is  not  possible. 

A  comparison  of  precipitation  forecasts  from  the  NGM  and  NORAPS  shows  that  the 
NGM  was  much  more  conservative  in  its  prediction  (not  sht^wn).  From  24  -  36  h,  the 
NGM  forecast  a  maxima  of  1"  (2.4  cm) .  while  Simulation  1  predicted  a  maxima  of  4.5" 
(10.9  cm)  and  Simulation  2  a  maxima  of  3.5”  (8.5  cm).  Precipitation  forecasts  for  the 
Spectral  model  were  not  available. 


VI.  MODEL  DIAGNOSTICS 


In  the  previous  section  it  was  observed  that  w'hile  the  48  h  prediction  by  the  NORAPS 
model  was  very  good  in  terms  of  central  pressure,  the  36  h  forecast  (valid  at  14/0000)  was 
16  mb  too  deep.  This  was  largely  the  result  of  intense  convection  early  in  the  simulation. 
The  model  run  with  no  precipitation  showed  only  a  small  error  in  central  pressure  at 
14/0000.  This  section  will  investigate  the  effects  of  convection  in  the  model  simulations. 
In  all  figures  of  model  data,  the  actual  sigma  level  values  were  used  to  construct  cross 
sections  and  Skew-T  diagrams. 

A.  TIME  SECTIONS 

Section  3  described  the  Kuo  .scheme  of  convection  parameterization  which  is  used  in 
the  NORAPS  model.  The  effective  result  of  convection  is  to  transport  latent  heat  and 
moisture  from  low  levels  (i.e.,  the  PBL)  upward.  But  before  convection  in  the  model 
begins,  four  criteria  must  be  met  (reviewcti  in  Section  3). 

It  was  already  shown  in  section  .s  that  the  convective  precipitation  was  occurring  in 
areas  with  convergence  at  1000  mb.  implying  thtit  the  criteria  for  low-level  moisture  con¬ 
vergence  had  been  met.  To  see  if  the  other  three  criteriti  were  correctly  satisfied,  time 
sections  w'ere  constructed  at  ERICA  drifting  buoy  loctitions  in  which  intense  convective 
precipitation  was  falling  in  the  model.  Figures  23  -  25  show'  the  equivalent  potential 
temperature,  relative  humidity,  and  mixing  ratio  time  sections  for  Simulation  1  at  33°N 
7(EN  (buoy  1 1369).  The  time  sections  clearly  show  th:it  the  Kuo  scheme  is  accomplishing 
its  goal.  The  profile  of  0^  indicates  a  deep  layer  (surface  to  850  mb)  of  conditional  insta¬ 
bility  due  to  cold  air  from  the  New  England  anticyclone  flowing  over  the  warm  Gulf 
Stream  w  aters.  The  air  is  also  nearly  satiir.'ited  up  to  650  mb  (Fig.  24).  Thus,  the  neces- 


sary  conditions  for  convection  outlined  in  section  3  appear  to  exist.  The  potential  tempera¬ 
ture  time  section  (not  shown)  indicates  weak  static  stability  below  850  mb  with  stronger 
stability  above  this  layer.  Thus,  the  conditional  instability  (d0g/dz  <  0)  in  Fig.  23  is  large¬ 
ly  due  to  the  moisture  profile. 

From  charts  of  convective  precipitation  (see  Figs.  14  and  15b),  it  was  determined  that 
very  light  precipitation  from  convection  began  in  NORAPS  between  12  and  18  h  into  the 
forecast  simulation  at  the  33°N  70°W  grid  point.  From  hours  18  to  24,  approximately  5 
cm  of  precipitation  fell  at  this  location.  Figures  23  and  25  indicate  that  the  convective 
parameterization  began  to  take  effect  after  about  the  12  h  mark,  since  the  upper  levels  (700 
-  500  mb)  begin  to  moisten  creating  a  conditionally  neutral  layer  aloft  as  shown  by  the 
vertical  0^  lines.  The  low-level  conditional  instability  also  begins  to  be  removed  after  12 
h.  A  sharp  jump  in  the  moisture  profile  occurs  18  h  into  the  simulation  (Fig.  25).  the  same 
time  at  which  the  significant  convective  precipitation  begins.  The  entire  column  from  the 
surface  to  500  mb  is  nearly  neutral  by  this  time.  The  low-level  moisture  has  been  moved 
upwards  as  evidenced  by  the  pocket  of  greater  than  00%  humidity  in  Fig.  24  and  the  drier 
70%  contours  below.  After  hour  21.  the  upper  levels  begin  to  dry  rapidly  (Fig.  24)  indicat¬ 
ing  the  intrusion  of  drier  air  from  the  west,  while  the  Itnv  levels  begin  to  moisten  as  evi¬ 
denced  by  the  90%  contour,  setting  up  another  conditionally  unstable  situation.  The  cause 
of  this  moistening  is  two-fold:  evaporation  of  the  precipitation  into  the  drier  low-level  air 
and  a  shift  in  the  low'-level  w'ind  direction  from  northeasterly  to  sttutherly  (not  showm)  as 
the  warm  front  passes  by,  advecting  mtue  moist  air  from  the  south. 

Figure  26  is  a  north-south  cross  section  from  a  study  by  Reed  and  Albright  (1986)  of 
an  explosively  deepening  cyclone  in  the  eastern  North  Pacific.  Note  hinv  the  equivalent 
potential  temperature  contours  are  oriented  vertically  where  convection  is  occurring.  The 
0J.  lines  in  the  current  study  (Fig.  23)  are  also  oriented  vertically  after  18  h  as  significant 
convective  precipitation  began  to  fall.  In  the  cold  air  to  the  north  (left)  in  Fig.  26.  potential 
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instability  is  indicated  by  the  d0g/dz  <  0  up  to  700  mb.  This  is  similar  to  the  left  side  of 
Fig.  23  where  the  cold  air  from  the  New  England  anticyclone  is  moving  over  the  warm 
water. 

A  comparison  between  the  air  temperatures  measured  by  the  buoy  and  the  analysis 
temperatures  at  each  model  level  was  performed  (not  shown).  This  comparison  revealed 
that  the  lowest  model  level  (21  m)  was  1.4°C  warmer  than  the  buoy.  Gjnsidering  that  the 
New  England  anticyclone  was  advecting  cold  air  over  the  warm  water,  one  would  expect 
to  find  colder  temperatures  above  the  surface  layer.  This  anomaly  was  also  observed  at 
other  buoy  locations  in  the  same  region.  The  temperature  difference  became  larger  as  the 
simulation  progres.sed.  The  cause  of  this  warming  is  mo.st  likely  the  warm  ocean  waters. 
But  a  comparison  of  the  measured  sea-surface  temperatures  (SST's)  and  the  model  SST 
analysis  showed  good  agreement.  This  appears  tci  imply  that  the  NORAPS  model  might 
possibly  be  modifying  the  cold  air  too  rapidly.  This  warming  by  the  warm  ocean  destabi¬ 
lizes  the  atmosphere,  as  seen  in  Fig.  23. 

Similar  time  sections  were  constructed  for  Simulations  2-4.  It  was  pointed  out  in 
section  5  that  Simulations  2  and  3  did  not  differ  significantly  from  Simulation  1  in  terms  of 
the  .sea  level  pressure  forecast,  but  that  the  convective  and  stable  precipitation  predictions 
were  different.  The  precipitation  charts  indicate  that  the  total  rainfall  at  the  diagnostic  grid 
point  for  Simulation  2  (KU074)  was  much  le.ss  than  Simulation  1  (KU065),  («  1.5  cm 
compared  with  «  5.0  cm).  The  KU074  scheme  intentionally  decreases  the  moistening 
when  saturation  is  nearly  reached,  which  decreases  the  timount  of  precipitation.  But  as  the 
moistening  is  decreased,  the  scheme  c('rrespondingly  increases  the  sensible  heating. 
Additionally,  stable  precipitation  is  increased  in  the  KU074  simulation.  Thus,  the  overall 
heating  profile  is  nearly  identical  to  Simuhition  1  (KU065). 

This  difference  between  the  two  convective  parameterization  schemes  is  shown  in  the 
time  sections.  Figures  27  and  2S  are  the  equivalent  potential  temperature  and  mixing  ratio 


time  sections  for  Simulation  2.  Comparing  these  to  Figs.  23  and  25,  both  simulations  look 
nearly  identical  for  the  first  15  h.  But  the  conditional  instability  is  removed  more  quickly 
in  Simulation  1  as  shown  by  comparing  the  vertical  0^  contours  in  Figs.  23  and  27.  The 
moistening  process  is  also  slower  in  the  second  simulation,  which  is  noted  by  comparing 
the  8  g/kg  contour  in  Figs.  25  and  28.  The  first  simulation  (KU065)  moistens  the  column 
more  quickly  than  Simulation  2  (KU075),  which  removes  the  conditionally  unstable  layer 
faster  and  causes  more  precipitation.  Thus,  the  intense  convection  in  Simulation  2  is 
delayed  as  compared  to  Simulation  1. 

As  would  be  expected,  the  time  sections  for  Simulation  4  (no  precipitation)  are  con¬ 
siderably  different  from  the  other  three  simulations.  Figures  29  and  30  show  the  equiva¬ 
lent  potential  temperature  and  the  mixing  ratit'  respectively.  The  initial  iow-ltv’el  condi¬ 
tional  instability  exists  as  in  Simulations  1  and  2.  However,  with  the  convective  parame¬ 
terization  turned  off,  the  model  has  no  way  to  relieve  this  instability.  Instead,  the  condi¬ 
tionally  unstable  layer  deepens.  Alstg  note  the  marked  differences  in  the  moisture  profiles 
of  the  twii  simulations,  with  Simulation  4  indicating  only  ;i  small  amount  of  moistening  in 
the  low  levels,  and  none  aloft.  During  this  period  { 15  to  27  h).  the  i(^w-level  winds  shift 
from  northeasterly  to  easterly,  again  suggesting  that  ads  ection  w  as  the  primary  cause  of 
the  moistening  observed  in  this  simulation. 

B.  IMPORTANCE  OF  INITIAL  CONDITIONS 

From  a  qualitative  analysis,  the  ccmvective  parameterization  scheme  appears  to  be 
performing  correctly.  The  model  atmosphere  possesses  the  necessary  eruditions  to  cause 
convection  and  the  Kuo  scheme  relieves  the  pr'tential  instability  by  transpc'rting  low  -level 
latent  heat  and  moisture  upward.  However,  the  high  values  of  convective  precipitation 
noted  in  section  4  indicate  that  the  model  is  perhrrming  incorrectly  in  its  prediction  of 


convection.  Since  the  predicted  temperature  and  moisture  play  a  large  role  in  the  cumulus 
parameterization,  these  fields  are  further  investigated. 

The  time  sections  in  section  5a  were  coincident  with  a  drifting  buoy  in  the  location  of 
predicted  convective  precipitation.  However,  this  did  not  allow  for  comparison  with  actual 
data  at  levels  above  the  surface.  In  order  to  achieve  this,  a  cross  section  was  constructed 
using  four  Air  Force  dropsondes  taken  at  13/0000.  Figure  31  is  the  equivalent  potential 
temperature  and  mixing  ratio  cross  sections  along  the  line  indicated  on  Fig.  14.  The  low- 
level  potential  instability  created  by  the  advection  of  cold  air  over  warm  water  is  clearly 
indicated  by  the  "tongue"  of  low  0^  values  near  900  mb.  The  corresponding  NORAPS  12 
h  forecast  cross  section  for  Simulation  1  is  shown  in  Fig.  32.  NORAPS  also  possesses  a 
layer  of  potential  instability.  From  Fig.  14  it  is  noted  that  convective  precipitation  has 
already  been  occurring  in  the  region  of  the  cross  section.  Thus,  some  of  the  instability  has 
already  been  relieved,  especially  in  the  region  between  dropsonde  3  and  4.  A  similar  cross 
section  6  h  earlier  (not  shown)  does  in  fact  show  that  the  model  atmosphere  was  more 
unstable. 

Comparison  of  Figs.  31  and  32  reveals  a  disturbing  fact:  the  model  atmosphere  is 
much  moister  in  the  low  levels,  typically  by  2  g/Tcg.  Analysis  and  6  h  cross  sections  (not 
shown)  show  little  change  in  the  moisture  field.  Thus,  it  appears  NORAPS  was  initialized 
with  a  poor  moisture  field.  As  was  previously  mentioned.  NORAPS  does  not  perform  an 
actual  moisture  analysis.  Instead,  it  uses  its  previous  12  h  moisture  forecast  as  the  analy¬ 
sis.  Obviously,  this  could  result  in  a  situation  such  as  that  in  the  present  study. 

Figure  33  shows  a  comparison  of  the  actual  Air  Force  dropsonde  3  and  the  corre¬ 
sponding  NORAPS  12  h  forecast  sounding  (denoted  by  /-\F03  in  Figs.  31  and  32).  While 
the  overall  soundings  look  similar,  closer  inspection  reveals  that  NORAPS  is  both  moister 
and  warmer  than  the  observed  atmosphere.  Also  note  that  the  observed  inversion  between 
850  -  700  mb  (Fig.  33a)  is  not  predicieti  by  the  model  (Fig.  33b).  The  result  is  a  more 
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unstable  sounding  in  the  model,  as  indicated  by  the  stability  indices.  Similar  soundings  for 
all  other  (18)  dropsondes  show  that  the  model  is  more  unstable  than  the  real  atmosphere. 
Tills  is  somewhat  surprising  since  the  actual  atmosphere  "looked"  more  unstable  than  the 
model  in  the  0g  cross  sections  (Pigs.  31a  and  32a).  ALso,  convective  activity  in  this  region 
had  already  served  to  stabilize  the  model  atmosphere  somewhat.  Clearly,  the  model  is  not 
doing  a  good  job  of  predicting  the  temperature  and  moisture  fields  in  the  lower  levels. 
And  since  the  analysis  fields  (not  shown)  are  even  more  unstable,  it  is  suspected  that  the 
initial  NORAPS  fields  poorly  represented  the  actual  atmosphere. 

In  order  to  determine  if  the  poor  representation  of  the  temperature  and  moi.sture  was 
due  the  lack  of  data,  comparisons  were  made  between  the  model  and  the  12/1200  sound¬ 
ings  (not  shown)  at  West  Palm  Beach,  FL  (72203).  Charleston,  SC  (72208),  Cape  Haiteras. 
NC  (72304),  Nassau,  Bahamas  (78073).  and  Bermuda  Island  (78016).  In  all  cases  but 
Bermuda,  the  model  analysis  was  nearly  identical  to  the  observaition  data.  As  with  the  Air 
Force  dropsondes.  NORAPS  w'as  too  vxiim  and  moist  ;it  Bermuda.  Therefore  it  appears 
that  the  NORAPS  analysis  gives  a  correct  depictitm  of  the  titmosphere  when  data  are  avail¬ 
able,  but  may  be  poor  over  the  data  sparse  ocean  regions.  Model  12  h  forecast  errors  are 
passed  on  to  the  next  analysis  through  the  first  guess  field  if  there  is  no  data  to  refute  it. 
Thus,  a  poor  12  h  forecast  can  result  in  a  ptntr  analysis. 

It  was  mentioned  in  section  3  that  for  the  present  study,  four  12  h  forecasts  had  been 
made  in  order  to  simulate  the  update  eycie  m  an  ('peratii'ital  enx  ironment.  Hodur  (1987) 
showed  that  the  inclusion  of  a  12  h  update  cycle  improved  NORAPS  forecasts.  Thus,  the 
NORAPS  analysis  could  differ  from  the  gh'bal  mcnlel  (NOO.APS)  analysis  because  of  the 
use  of  the  previous  12  h  forecast  :is  the  first  guess  field.  Also  in  section  3  it  was  pointed 
out  that  neither  NORAPS  nor  NOOAPS  antdyze  moisture.  Rather,  the  initial  moisture 
fie.d  for  a  model  simulation  is  the  previous  12  h  lorecast.  Surface  tempertitures  measured 
by  ships  and  buens  are  also  not  used  in  the  tinalyMs. 


To  investigate  the  impact  of  these  procedures.  Fig.  34  compares  the  2  m  air  tempera¬ 
ture  analysis  from  NORAPS  (34a)  and  NOGAPS  (34b)  with  actual  ship  and  buoy  observa¬ 
tions.  The  ERICA  buoy  data  shown  in  this  figure  was  not  used  in  the  model  analysis. 
Note  that  south  of  35°N,  NORAPS  (Fig.  34a)  is  generally  warmer  than  the  observations, 
by  as  much  as  6°C.  The  NOGAPS  analysis  (Fig.  34b)  shows  much  less  detail  in  the 
temperature  field,  especially  near  the  coastline.  However,  the  analysis  appears  to  be  in 
much  better  agreement  with  the  observed  temperatures,  although  the  ERICA  buoys  near 
32°N  are  still  colder  than  the  analysis.  Since  the  temperature  is  not  analyzed  directly  by 
either  model,  these  differences  are  largely  a  result  of  the  previous  12  h  forecast. 

Figure  35  is  the  initial  1000  mb  vapor  pressure  (moisture)  fields  for  the  NORAPS 
(35a)  and  NOGAPS  (35h)  models.  The  NORAPS  model  again  shows  more  detail  in  the 
1000  mb  moisture  field  than  does  the  global  model.  However,  the  regional  model  is 
considerably  moister,  especially  south  of  35°N  (note  the  Icication  of  the  15  mb  contour  as  it 
crosses  70°W).  Comparison  of  actual  moisture  values  at  12/1200  (not  shown)  to  the 
NORAPS  and  NOG.APS  analyses  show  that  NOGAPS  is  more  much  more  realistic,  and 
may  in  fact  be  slightly  too  dry.  The  better  moisture  analyses  in  NOG.APS  are  actually 
better  12  h  forecasts,  since  there  is  no  actual  analysis  of  moisture. 

The  causes  for  the  differences  between  the  NORAPS  ;ind  NOGAPS  low-level  temper¬ 
ature  and  moisture  fields  are  beyond  the  scope  of  this  study.  As  mentioned  in  section  6a. 
one  possibility  is  that  as  the  continental  pohir  air  mt'ves  over  the  Gulf  Stream  waters,  the 
NORAPS  model  is  modifying  the  air  too  rapidiv,  thus  mtiking  it  warmer  and  moister  in  the 
low  levels.  The  drier  NOGAPS  fields  might  :dso  be  the  result  of  different  model  physics, 
such  as  the  use  of  the  Arakawa-Sehiibert  ( 1M74)  scheme  bm  convection  (Rosmond'^. 

4.  Dr.  Thomas  H.  Rosmond  is  the  head  of  the  Rredietiim  Systems  bnineh  of  the  Ntival  Oeetinogniphic  & 
Atmospheric  Research  I.abor;itory,  Atmospheric  Directorate.  Monterey.  ('A. 


personal  communication),  rather  than  the  Kuo  (1974)  scheme  used  in  the  NORAPS  model 
for  the  current  study. 

C.  SIMULATION  WITH  NOGAPS  INITIAL  FIELDS 

In  order  to  test  the  hypothesis  that  the  excess  temperature  and  moisture  in  the 
NORAPS  initial  fields  are  adversely  effecting  the  model  simulations,  a  fifth  simulation 
was  carried  out  where  the  model  was  initialized  with  the  global  analysis  fields  (i.e.,  no 
regional  update  cycle).  The  model  for  the  this  simulation  used  the  KU074  convection 
scheme.  Additionally,  minor  improvements  were  made  to  the  boundary  conditions,  but 
these  have  relatively  small  influence  on  the  forecast  (Liou,  personal  communication). 

A  subjective  comparison  of  the  NOGAPS  surface  and  500  mb  analyses  with  the  fields 
used  to  start  the  first  four  simulations  showed  scmie  minor  differences,  but  these  were 
primarily  out  of  the  region  of  interest.  The  important  weak  disturbance  over  Florida  was 
depicted  in  both  analyses  by  the  8  s'*  vcuticity  contour.  Thus,  the  major  difference  be¬ 
tween  this  simulation  and  Simulation  2  is  the  ini'ial  low-level  temperature  and  moisture 
fields. 

The  cumulative  convective  precipitation  from  12  to  18  h  is  shown  in  Fig.  36.  Com¬ 
paring  this  with  Fig.  19  from  Simulation  2.  it  is  seen  that  the  drier  initial  moisture  iicid  has 
had  an  effect  on  the  convective  precipitatic'n.  both  in  amount  and  distribution.  The  corre¬ 
sponding  sea  level  pressure  fields  (not  shown)  differ  only  by  ]  mb.  Likewise,  the  500  mb 
fields  show  little  difference,  although  the  new'  simuhition  moves  the  .short  wave  somew'hat 
faster.  Still,  the  model  is  developing  the  w  rong  inititil  low.  But  by  hour  30  (Fig.  37),  an 
improvement  is  noted.  Since  the  inititil  low  at  this  time  is  not  as  deep  as  in  Simulation  2 
(998  mb  vs  991  mb  in  Simulation  2,  not  show  n).  NORAPS  attempts  to  develop  the  sec¬ 
ondary  low  to  the  northwest.  The  subjective  analysis  from  Chalfant  (1989)  at  this  time 
(not  shown)  shows  the  ob.served  secimd  low  in  gocn!  agreement  with  the  forecast  position 
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and  central  pressure.  However,  the  primary  low  has  a  large  error  in  both  position  and 
central  pressure,  although  it  is  closer  to  the  ob.served  1002  mb  value  than  Simulation  2. 

By  14/0000,  the  two  simulations  are  considerably  different.  Figure  38  is  the  surface 
pressure  forecast  for  Simulation  5.  (The  forecast  by  Simulation  2  at  this  time  is  identical 
to  Simulation  1  (Fig.  16)  and  is  thus  not  shown).  The  989  mb  low  is  still  in  error  by  7-9 
mb.  But  this  is  a  considerable  improvement  over  Simulations  1-3.  Note  that  the  low  in 
Simulation  5  is  moving  much  more  rapidly.  Anthes  et  al.  (1983)  point  out  that  latent 
heating  reduces  the  cyclone's  speed  of  movement.  As  Table  2  indicates,  this  simulation 
produced  much  less  precipitation,  thus  agreeing  with  this  finding. 

Table  2  also  reveals  the  effect  of  the  reduced  precipitation.  The  48  h  central  pressure 
forecast  is  971  mb,  in  close  agreement  with  the  obsers  ed  968  mb.  Note  tiiat  it  appears  that 
the  amount  of  predicted  deepening  is  related  to  the  total  amount  of  precipitation,  not  the 
convective  or  stable  amounts  alone.  KU074  reduced  the  amount  of  convection  from  that 
of  KU065,  but  the  stable  precipitation  incretised  resulting  in  :i  similar  amount  of  latent 
heating.  Only  when  both  types  td'  precipitation  were  reduced  was  the  deepening  rate 
improved  (see  Fig.  21). 

Therefore,  it  appears  that  the  initial  h'w-level  temperature  and  moisture  fields  play  a 
significant  role  in  the  de\  elopment  of  the  cyclone.  The  cumulus  parameterization  is  sensi¬ 
tive  to  small  changes  in  the  moisture  field  ;md  correspondingly  varies  the  amount  of  pre¬ 
cipitation,  and  thus  the  timount  of  latent  heat  release.  The  result  is  a  difference  in  both  the 
track  and  intensity  of  the  predicteil  cyclone,  with  both  improved  over  the  earlier  simula¬ 
tions  w'ith  the  NORAPS  update  cycle.  Phe  ex.act  ctiuse  of  the  excessive  heat  and  moisture 
in  the  NORAPS  analyses  is  beyond  the  sco|ic  (d  the  current  study,  but  is  most  likely  linked 
to  the  boundary  layer  parameterization  (mixing)  and  vertical  grid  resolution  (see  Table  1). 
Again  it  should  be  noted  that  this  is  just  one  case  :md  th:it  similar  studies  mtty  show  results 
which  support  or  contradict  these  findings. 


VII.  SUMMARY  AND  CONCLUSIONS 


Considerable  work  has  been  done  in  the  area  of  numerical  prediction  of  rapid  oceanic 
cyclogenesis.  Several  studies  have  focused  on  the  latent  heating  aspect  as  the  key  mecha¬ 
nism  controlling  the  deepening  rate.  Yet,  in  nearly  all  of  these  studies,  the  deepening  rate 
was  typically  less  than  observed.  While  the  operational  NORAPS  48  h  forecast  valid  at 
14/1200  UTC  appeared  quite  good,  a  closer  in.spection  reveals  that  the  intermediate  fore¬ 
casts  were  rather  poor,  with  NORAPS  over-forecasting  the  development  of  the  ERICA 
lOP-2  cyclone. 

An  investigation  into  this  numerical  forecast  showed  that  the  NORAPS  model  incor¬ 
rectly  developed  a  low  too  early  in  the  forecast  due  to  pr(4'»lems  in  the  convective  parame¬ 
terization.  Furthermore,  the  deveh'pment  of  this  low  was  uh-*  rapid,  with  an  over-forecast 
central  pressure  of  18  mb  at  14/0000  UTC.  Unrealistically  large  amounts  of  convective 
precipitation  appear  to  be  responsible  for  the  overdevelopment,  as  the  convection  generat¬ 
ed  spurious  amounts  of  vorticity  at  500  mb  with  no  appreciable  change  in  the  height  pat¬ 
tern. 

Varying  the  cumulus  parameterization  scheme  decreased  the  amount  of  convective 
precipitation,  but  since  the  stable  precipitation  c('rrcspimdingly  increased,  tb.is  resulted  in 
very  little  change  in  the  surface  pressure  :md  500  mb  height  forecasts.  Danard  (1985) 
found  a  direct  correlation  between  the  mtiximum  amount  of  convective  precipitation  and 
the  predicted  deepening.  This  fact  was  not  observed  in  the  present  study.  Kuo  and  Low- 
Nam  (1990)  emphasized  the  importance  <4' warm  frontal  sttible  precipitation  to  the  rate  of 
deepening.  They  found  that  the  simulations  using  the  Arakawa-Schubert  convective 
parameterization  were  superior  to  those  using  Kuo.  They  concluded  that  this  was  primari¬ 
ly  the  result  of  the  former  scheme  precipitating  much  less  in  convection  than  the  latter. 


leaving  more  available  moisture  for  the  stable  precipitation.  Stable  precipitation  is  typical¬ 
ly  more  efficient  than  convective  precipitation  at  deepening  a  cyclone  since  it  occurs  lower 
in  the  atmosphere.  As  Anthes  et  al.  (1983)  points  out,  the  lower  the  level  of  maximum 
heating,  the  more  intense  the  resulting  surface  cyclone.  For  the  current  study,  the  total 
amount  of  precipitation  was  the  only  link  to  the  predicted  deepening  rate.  However,  the 
importance  of  the  location  of  the  precipitation,  both  horizontal  and  vertical,  cannot  be 
overlooked. 

In  order  to  verify  that  the  latent  heat  release  was  the  cause  of  the  overforecast  deepen¬ 
ing,  a  simulation  with  no  precipitation  was  run.  This  alteraticm  dramatically  slowed  the 
development  of  the  initial  cyclone,  but  tailed  to  adeqiiatelv  predict  the  explosive  deepening 
w'hich  was  observed.  The  underprediction  of  the  low  prior  t(^  the  rapid  deepening  phase 
indicates  that  precipitation  w’as  important  in  the  earlv  sttiges. 

Detailed  diagnostics  were  performed  to  tmalyze  the  role  of  the  convective  precipita¬ 
tion.  It  w'as  determined  that  the  cumulus  ptirameteriztition  scheme  was  working  properly. 
How'ever.  comparisons  of  model  antilyses  .and  forectists  to  ob.served  data  indicate  that  the 
model  was  poorly  initialized,  especitilly  in  low  levels  tner  the  ocetin.  resulting  in  a  great 
deal  more  potential  instability.  The  comptirisims  w  ith  HRICA  dropsimdes  also  showed 
that  the  model  forecasts  w  ere  too  w  arm  in  the  low  levels,  w  hich  mav  be  a  result  of  a  poor 
surface  flux  parameterization.  It  w  as  reci'gnized  that  the  lack  of  a  moisture  anaUsis  also 
might  be  partly  responsible  for  the  poor  initialization.  In  ('rder  to  investigate  this,  a  fifth 
simulation  was  conducted  which  used  the  global  anaivses  as  imtitil  fields.  These  global 
fields  were  found  to  be  drier  and  more  realistic.  The  resulting  forecasts  were  improved, 
both  in  sea  level  pressure  and  precipitation  amounts.  altlKuigh  the  model  still  deveh^ped  the 
wrong  low',  and  then  failed  to  deepen  it  enough. 

A  comparison  with  tw'o  NMC  models,  the  Nested  Clrid  Model  and  the  Global  Spectral 
Model,  showed  simihir  deficiencies  t<'  those  found  in  NORAPS.  Both  miulels  developed 


the  initial  low  and  tailed  to  resolve  the  secondary  development  off  the  Virginia  -  North 
Carolina  coastline.  While  the  NGM  model  was  l(i  mb  too  shallow  at  14;  1200,  the  spectral 
model  actually  overforecast  the  development  by  4  mb,  similar  to  the  operational  NORAPS 
forecast.  Both  models  did  very  well  in  predicting  the  48  h  location  of  the  cyclone. 

This  study  agrees  with  previous  findings  on  the  importance  of  latent  heating  in  the 
development  of  wintertime  oceanic  extralropical  cyclones.  The  quality  of  numerical 
forecasts  relies  on  the  quality  of  the  initial  data.  In  data  sparse  regions,  model  forecast 
errors  often  propagate  through  successive  forecasts  via  the  first  guess  field  unless  there  are 
available  data  to  correct  the  analysis.  Parameterizations  f(U  convection  and  surface  fluxes 
are  very  sensitive  to  details  in  the  miHlel  atmosphere.  Incorrect  specification  of  this  initial 
atmosphere  can  have  an  ad\  erse  effect  on  the  forecasts  even  if  the  parameterization 
schemes  are  correct.  Thus,  the  common  lack  (d'data  is  :m  important  problem  in  the  correct 
prediction  of  rapidly  developing  cycUmes. 

Since  the  initial  moisture  field  tippears  to  h.'ive  had  ;i  large  effect  on  the  NOR.APS 
forecasts,  future  work  should  be  done  in  this  ;»rea.  The  incc^rpirration  of  an  actual  moisture 
analysis,  rather  than  the  pres  ii'us  12  h  forecast,  cmikl  have  :t  significant  impact  on  explo¬ 
sive  cyclogenesis  forecasting.  Additiontilly.  it  was  noted  that  the  low  -level  temperatures  in 
NORAPS  were  often  too  warm  in  data  sparse  regions,  further  adriing  tir  the  potential  insta¬ 
bility  in  the  model.  The  current  studv  used  onl)  tin'se  ('bsers ations  which  were  ('per.ition- 
ally  available  at  the  time  (4  lOP-2.  The  extra  buox,  raw  insonde.  and  dropsonde  data  col¬ 
lected  as  a  part  of  ERICA  should  be  mtegr.iterl  into  the  NOR.APS  analyses  to  iinestigate 
the  importance  c.f  these  data  in  numerical  forecasting. 
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Fig.  1.  FNOC  500  mb  analysis  for  1200  UTC  12  December  19S8.  Geopotential  height 
(solid)  IS  in  dekameters  (+64  represents  564  dam)  and  temperature  (dashed)  is  in  '"C. 


f-'ig.  3.  f'NOC  surtacc  analysis  for  OOOO  ll'I'C'  13  December  19NiS.  Sea  le\el  pressure 
contours  are  in  mb  minus  lOOO  nrb  (  +  20  represents  1020  mb).  Surface  re|)orl  iiuls  aie 
in  kt  and  temperatures  are  in  °C.  See  text  for  explanation  of  circled  buoy. 
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As  ill  Fig,  3  except  tor  1200  UTC'  13  December  1988. 
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in  Fig.  1  except  for  1200  UTC  13  December  1988. 
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Fig.  6.  Subjective  surface  analysis  for  0000  UTC  14  December  1988.  Surface  pres.sure 
contour  interval  is  2  mb.  Surf  ace  observations  are  of  sea  level  pressure  (125  =  1012.5  mb), 
temperature  (in  °C)  and  winds  (in  kt).  From  Chalfant  (1989). 
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Fig.  8.  As  in  Fig.  6  except  for  1200  UTC  14  December  1988. 


Fig.  9.  As  in  Fig.  1  except  for  1200  UTC  14  December  1988. 
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Fig.  11.  12  h  forecast  of  sea  level  pressure  (contour  interval  4  mb)  and 
1000-500  mb  thickness  (dashed,  contour  interval  60  m)  for  Simulation  1. 
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Fig.  12a.  Initial  1000  mb  streamlines  for  Simulation  1. 
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Fig.  12b.  Initial  500  mb  geopotential  heights  (contour  interval  60  m)  and 
absolute  vorticity  (dashed  lines,  contour  interval  2  s'^)  for  Simulation  1. 


Fig.  13b.  As  in  Fig.  12b  except  tor  12  h  torecast. 
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Fig.  14.  12  h  forecast  of  cumulative  12  h  convective  precipitation  (contour 
interval  1  cm).  Solid  line  indicates  cross  section  location  in  Figs.  31  &  32. 
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Fig.  15a.  18  h  forecast  of  500  mb  geopoteniial  heights  (contour  interval  60 
m)  and  absolute  vorticity  (dashed  fines,  contour  interval  2  s'^)  for  Simula- 


Fig.  15b.  18  h  forecast  of  cumulative  6  h  convective  precipitation  (contour 
interval  1  cm)  for  Simulation  1. 
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Fig.  16.  As  in  Fig.  1 1  except  for  36  h  forecast. 
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Fig.  17b.  As  in  Fig.  12b  except  for  48  h  forecast  for  Simulation  1. 
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Fig.  18.  As  in  Fig.  15b  except  for  Simulation  2. 
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Fig.  19a.  As  in  Fig.  13a  except  for  Simulation  4. 
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Fig.  19b.  As  in  Fig.  13b  except  for  Simulation 


Central  Preseure  vs  Ttee 
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Fig.  21.  Observed  and  f  jrecast  central  pressures  (CP)  in  mb.  Observed  CP  taken  trom 
Chalfant  (1989).  Forecasts  are  for  Simulations  1-5  (see  text  for  explanation  on  simuki- 
tions). 
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Fig.  22a.  NMC  Global  Spectral  Model  36  h  forecast  of  sea  level  pressure 
(contour  interval  4  mb)  and  1000-500  mb  thickness  (dashed  lines,  contour 
interval  60  m). 
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Fig.  22b.  NMC  Global  Spectral  Model  36  h  forecast  of  500  heighi.'<  (con¬ 
tour  interval  60  m)  and  vorlicity  (dashed  lines,  contour  interval  2*  ). 
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Fig.  28.  As  in  Fig.  25  except  for  Simulation 


Fig.  31.  Vertical  cross  section  of  (a)  equivalent  potential  temperature  (contour  interval  2 
K)  and  (b)  mixing  ratio  (contour  interval  1  g/kg)  of  four  Air  Force  dropsondes  at  UtKK) 
UTC  13  December  1988.  See  Fig.  12  for  location  of  cross  section. 
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Fig.  32.  As  in  Fig.  31  except  for  Simulation  1  12  h  forecast. 
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Fig.  33.  Skew-T  Log  P  diagram  of  (a)  Air  Force  03  dropsonde  and  (b)  corre¬ 
sponding  NORAPS  12  h  forecast  valid  at  0000  UTC  13  December  1988. 
Temperatures  are  in  °C.  Stability  indices  at  right  are  Showalier  Index  (SI). 
Lifted  Index  (LI),  K  Index  (KI)  and  Totals  Index  (TTl). 


Fig.  34a.  1200  UTC  12  December  1988  analysis  of  2  m  air  temperature  by 
NORAPS.  Differences  between  observed  ship  and  buoy  temperatures  and 
the  model  analyses  are  also  plotted.  A  positive  value  indicates  that  the 
analysis  is  w'armer  than  the  observation.  Underlined  numbers  indicate 
ERICA  buoys. 


Fig.  34b.  As  in  Fig.  34a  except  for  NOGAPS  analysis. 


Fig.  35a.  1200  UTC  12  December  1988  analysis  of  1000  mb  vapor  pres¬ 
sure  (contour  interval  5  mb)  by  NORAPS. 
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